Avalanches and the Renormalization Group for Pinned Charge-Density Waves

The critical behavior of charge-density waves (CDWs) in the pinned phase is studied for applied fields increasing toward the threshold field, using recently developed renormalization group techniques and simulations of automaton models. Despite the existence of many metastable states in the pinned state of the CDW, the renormalization group treatment can be used successfully to find the divergences in the polarization and the correlation length, and, to first order in an $\epsilon = 4-d$ expansion, the diverging time scale. The automaton models studied are a charge-density wave model and a ``sandpile'' model with periodic boundary conditions; these models are found to have the same critical behavior, associated with diverging avalanche sizes. The numerical results for the polarization and the diverging length and time scales in dimensions $d=2,3$ are in agreement with the analytical treatment. These results clarify the connections between the behaviour above and below threshold: the characteristic correlation lengths on both sides of the transition diverge with different exponents. The scaling of the distribution of avalanches on the approach to threshold is found to be different for automaton and continuous-variable models.Comment: 29 pages, 11 postscript figures included, REVTEX v3.0 (dvi and PS files also available by anonymous ftp from external.nj.nec.com in directory /pub/alan/cdwfigs

Connectivity strategies to enhance the capacity of weight-bearing networks

The connectivity properties of a weight-bearing network are exploited to enhance it's capacity. We study a 2-d network of sites where the weight-bearing capacity of a given site depends on the capacities of the sites connected to it in the layers above. The network consists of clusters viz. a set of sites connected with each other with the largest such collection of sites being denoted as the maximal cluster. New connections are made between sites in successive layers using two distinct strategies. The key element of our strategies consists of adding as many disjoint clusters as possible to the sites on the trunk $T$ of the maximal cluster. The new networks can bear much higher weights than the original networks and have much lower failure rates. The first strategy leads to a greater enhancement of stability whereas the second leads to a greater enhancement of capacity compared to the original networks. The original network used here is a typical example of the branching hierarchical class. However the application of strategies similar to ours can yield useful results in other types of networks as well.Comment: 17 pages, 3 EPS files, 5 PS files, Phys. Rev. E (to appear

Vortex Dynamics and Defects in Simulated Flux Flow

We present the results of molecular dynamic simulations of a two-dimensional vortex array driven by a uniform current through random pinning centers at zero temperature. We identify two types of flow of the driven array near the depinning threshold. For weak disorder the flux array contains few dislocation and moves via correlated displacements of patches of vortices in a {\it crinkle} motion. As the disorder strength increases, we observe a crossover to a spatially inhomogeneous regime of {\it plastic} flow, with a very defective vortex array and a channel-like structure of the flowing regions. The two regimes are characterized by qualitatively different spatial distribution of vortex velocities. In the crinkle regime the distribution of vortex velocities near threshold has a single maximum that shifts to larger velocities as the driving force is increased. In the plastic regime the distribution of vortex velocities near threshold has a clear bimodal structure that persists upon time-averaging the individual velocities. The bimodal structure of the velocity distribution reflects the coexistence of pinned and flowing regions and is proposed as a quantitative signature of plastic flow.Comment: 12 pages, 13 embedded PostScript figure

Exact Phase Diagram of a model with Aggregation and Chipping

We revisit a simple lattice model of aggregation in which masses diffuse and coalesce upon contact with rate 1 and every nonzero mass chips off a single unit of mass to a randomly chosen neighbour with rate $w$. The dynamics conserves the average mass density $\rho$ and in the stationary state the system undergoes a nonequilibrium phase transition in the $(\rho-w)$ plane across a critical line $\rho_c(w)$. In this paper, we show analytically that in arbitrary spatial dimensions, $\rho_c(w) = \sqrt{w+1}-1$ exactly and hence, remarkably, independent of dimension. We also provide direct and indirect numerical evidence that strongly suggest that the mean field asymptotic answer for the single site mass distribution function and the associated critical exponents are super-universal, i.e., independent of dimension.Comment: 11 pages, RevTex, 3 figure

Order in extremal trajectories

Given a chaotic dynamical system and a time interval in which some quantity takes an unusually large average value, what can we say of the trajectory that yields this deviation? As an example, we study the trajectories of the archetypical chaotic system, the baker's map. We show that, out of all irregular trajectories, a large-deviation requirement selects (isolated) orbits that are periodic or quasiperiodic. We discuss what the relevance of this calculation may be for dynamical systems and for glasses

Evolution of Non-Equilibrium Profile in Adsorbate Layer under Compressive Strain

We investigate the time evolution of an initial step profile separating a bare substrate region from the rest of the compressively strained adsorbate layer near a commensurate to incommensurate transition. The rate of profile evolution as a function of the mismatch, coverage and the strength of the substrate potential are determined by Brownian molecular dynamics simulations. We find that the results are qualitatively similar to those observed for the Pb/Si(111) system. The anomalously fast time evolution and sharpness of the non-equilibrium profile can be understood through the domain wall creation at the boundary and its subsequent diffusion into the interior of the adsorbate layer.Comment: 6 pages, 7 figures, Tribology Letter

Vector lattice model for stresses in granular materials

A vector lattice model for stresses in granular materials is proposed. A two dimensional pile built by pouring from a point is constructed numerically according to this model. Remarkably, the pile violates the Mohr Coulomb stability criterion for granular matter, probably because of the inherent anisotropy of such poured piles. The numerical results are also compared to the earlier continuum FPA model and the (scalar) lattice $q$-model

Properties of layer-by-layer vector stochastic models of force fluctuations in granular materials

We attempt to describe the stress distributions of granular packings using lattice-based layer-by-layer stochastic models that satisfy the constraints of force and torque balance and non-tensile forces at each site. The inherent asymmetry in the layer-by-layer approach appears to lead to an asymmetric force distribution, in disagreement with both experiments and general symmetry considerations. The vertical force component probability distribution is robust and in agreement with predictions of the scalar q model while the distribution of horizontal force components is qualitatively different and depends on the details of implementation.Comment: 18 pages, 12 figures (with subfigures), 1 table. Uses revtex, epsfig,subfigure, and cite. Submitted to PRE. Plots have been bitmapped. High-resolution version is available. Email [email protected] or download from http://rainbow.uchicago.edu/~mbnguyen/research/vm.htm

Roughness at the depinning threshold for a long-range elastic string

In this paper, we compute the roughness exponent zeta of a long-range elastic string, at the depinning threshold, in a random medium with high precision, using a numerical method which exploits the analytic structure of the problem (`no-passing' theorem), but avoids direct simulation of the evolution equations. This roughness exponent has recently been studied by simulations, functional renormalization group calculations, and by experiments (fracture of solids, liquid meniscus in 4He). Our result zeta = 0.390 +/- 0.002 is significantly larger than what was stated in previous simulations, which were consistent with a one-loop renormalization group calculation. The data are furthermore incompatible with the experimental results for crack propagation in solids and for a 4He contact line on a rough substrate. This implies that the experiments cannot be described by pure harmonic long-range elasticity in the quasi-static limit.Comment: 4 pages, 3 figure

Behaviour of the energy gap near a commensurate-incommensurate transition in double layer quantum Hall systems at nu=1

The charged excitations in the system of the title are vortex-antivortex pairs in the spin-texture described in the theory by Yang et al which, in the commensurate phase, are bound together by a ``string''. It is shown that their excitation energy drops as the string lengthens as the parallel magnetic field approaches the critical value, then goes up again in the incommensurate phase. This produces a sharp downward cusp at the critical point. An alternative description based on the role of disorder in the tunnelling and which appears not to produce a minimum in the excitation energy is also discussed. It is suggested that a similar transition could also occur in compressible Fermi-liquid-like states.Comment: latex file, 17 page